Locally anisotropic momentum distributions of hadrons at freeze-out in relativistic 

heavy-ion collisions 
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A spheroidal anisotropic local momentum distribution is implemented in the statistical model of 
hadron production. We show that this form leads to exactly the same ratios of hadronic abundances 
as the equilibrium distributions, if the temperature is identified with a characteristic transverse- 
momentum scale. Moreover, to a very good approximation the transverse-momentum spectra of 
hadrons are the same for isotropic and anisotropic systems, provided the size of the system at 
freeze-out is appropriately adjusted. We further show that this invariance may be used to improve 
the agreement between the model and experimental HBT results. 
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I. INTRODUCTION 

Thermal analyses of hadronic yields in heavy-ion col- 
lisions [H-l25j have become a popular tool for the inter- 
pretation of the data collected at the AGS (ToUl^. SPS 
[THH, and RHIC energies (20^25). Most often, the suc- 
cessful results of such analyses are understood as the ev- 
idence for high-level thermalization of the hadronic sys- 
tems produced in such collisions. 

If the single- freeze-out scenario is assumed [22T - l24j . in 
addition to the hadron yields, the thermal approach may 
be used to describe other physical observables. In par- 
ticular, it has been successfully used to reproduce the 
RHIC transverse momentum spectra [26|, H3, collective 
flow [2f| [H, [23], femtoscopic observables |3fj| - l32| . and 
charge balance functions [H, [34[ . 

In the thermal approach, one usually obtains very good 
fits with just a few thermodynamic parameters, such as 
the temperature, T, the baryon and strangeness chemical 
potentials, /ib and /is, or in the extended approaches the 
quark fugacities, 7 S and 7 9 , and additional parameters 
describing the system's geometry and expansion (flow). 

Recently, the single-freeze-out model has been applied 
to the LHC data describing transverse-momentum spec- 
tra of hadrons produced in Pb+Pb collisions at the beam 
energy ^/snn = 2.76 TeV |35j. It has been shown that 
the model reproduces well the spectra of pions, kaons, 
and hyperons, however, it does not reproduce the pro- 
ton spectra. This results agrees with an earlier finding 
that the ratios of all hadron abundances measured at the 
LHC, except for protons, may be described within the 
thermal model 36] . 

In this paper we show that the agreement of the ther- 
mal models with the data does not necessarily imply the 
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fact that the system is locally equilibrated. We take 
into account possible differences in the local distribu- 
tion of transverse and longitudinal momenta (with the 
directions defined with respect to the beam axis). Us- 
ing the covariant Romatschke-Strickland form (RSF) for 
the anisotropic phase-space momentum distributions, we 
show that i) the ratios of hadron multiplicities are exactly 
the same as in the locally equilibrated system, if the tem- 
perature T is identified with a characteristic transverse- 
momentum scale A, ii) the transverse-momentum spectra 
of hadrons are changed, to a very good approximation, 
only by an overall factor which depends on the momen- 
tum anisotropy parameter appearing in RFS — since this 
change may be easily compensated by a change of the pa- 
rameters defining the size of the system, the transverse- 
momentum spectra may be also treated as insensitive to 
the discussed momentum anisotropy, finally, iii) the dis- 
cussed invariance of the spectra may be used to improve 
the model HBT results for pions. 

Generally speaking, our results demonstrate robust- 
ness of the thermal framework against specific variations 
of the model assumptions. This and other similar studies 
(see, for example, [37] ) show that the thermal approach 
is quite successful even in the situations where matter 
is out of equilibrium. Our results concerning the ratios 
of hadronic abundances agree qualitatively with the re- 
sult of Ref. [H|], where an approximate invariance of 
the ratios has been demonstrated with slightly different 
physical assumptions. 

We stress that we consider the anisotropy between lon- 
gitudinal and transverse momenta considered in the local 
rest frame of the fluid element. These two directions are 
defined always with respect to the beam axis. We do 
not consider the anisotropy of the transverse-momentum 
distributions in the (p x ,Py) plane, quantified by the har- 
monic flow coefficients v n . 

For simplicity, in the discussion of the spectra we con- 
sider boost-invariant and cylindrically symmetric models 
that have been implemented in the Monte-Carlo gener- 
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ator THERMINATOR [4J, |43|. We compare the results of 
the model calculations with the LHC data for Pb+Pb 
collisions at the beam energy ^/snn = 2.76 TeV. 



II. COOPER-FRYE FORMALISM AND 
COVARIANT ROMATSCHKE-STRICKLAND 
FORM 

In this paper we assume that primordial hadrons are 
emitted from the freeze-out hypersurface £ and their 
spectra may be calculated with the help of the Cooper- 
Frye formula 



dN 

E p — = / dV^x)p»f{x,p), 



(1) 



where E p is the particle's energy, cZE^ is the element 
of the freeze-out hypersurface, and f(x,p) is the phase- 
space distribution function of the emitted particles. For 
systems in local thermal equilibrium one assumes Fermi- 
Dirac (e = —1) or Bose-Einstein (e = +1) statistics and 
uses the distribution functions 



feq(x,p) = g { exp 



p-U - /i 



T 



(2) 



where g is the number of internal degrees of freedom, 
T is the freeze-out temperature, and /i is the chemical 
potential. Typically, the freeze-out conditions are defined 
by the fixed values of T and fi, hence these two quantities 
may be treated as constants in ([2}. 

The equilibrium distributions depend on the scalar 
product of the particle four-momentum p^ and the four- 
velocity of the fluid element. We use the standard 
parameterizations 

P 1 * = ( E P'P±>P\\) = (m±cosh.y,p ± ,m±smhy) (3) 
and 



U» = 7 (1, ""J 



(4) 



where 7 = (1 — v 2 ) -1 / 2 . 

In order to study the effects of the anisotropic distri- 
butions, we modify Eq. ([3} and use (39l - l4l1 | 



/rsf(^,p) = g { exp 



A 



(5) 

Here £ is the anisotropy parameter and A is a typical 
transverse-momentum scale characterizing the particles 
in the system. For local thermodynamic equilibrium 
£ = and A = T (from now on we set /i = 0). 

The four-vector appearing in ([5]) defines the direc- 
tion of the beam and has the structure 



V t = -y z {v z , 0,0,1), 7z = (l-« 2 ) 



'1/2 



(0) 



We note that the four-vectors and satisfy the 
normalization conditions 



U 2 = 1, V 2 = -l, U-V = 0. 



(7) 



In the local rest frame (LRF) of the fluid element, 
and have simple forms 

tf" = (1,0,0,0), ^ = (0,0,0,1), (8) 

and the distribution function ([5} is reduced to 

Jrsf(x,p) = g {exp 

(?) 

Neglecting the hadron mass (m = 0) and quantum statis- 
tics (e = 0) we recover the commonly used RSF (39j . 

III. RATIOS OF HADRONIC YIELDS 

Starting with the Cooper-Frye formula ([TJ we obtain 
the multiplicity of the hadron species i in the form 



N 4 



If the distribution function is taken as RSF, we have 
fi(x,p) — /rfs(p ' U 7 p ■ V), and the integral over the 
momentum in (|10[) yields the particle current 



d 3 p 
~E r , 



p^f^sip ■U,p-V) = n i (A, £(x)) U». (11) 



A few comments are in order now: i) we assume now that 
the freeze-out is defined by the fixed value of the hard 
scale A, hence, A is kept constant in (TTTl) . ii) we allow 
for the spacetime dependence of the anisotropy parame- 
ter £ but we assume that £ is the same for all hadronic 
species, hi) there is no term proportional to on the 
right hand of Eq. (|11[) due to quadratic dependence of 
the distribution function on p ■ V. 

The calculation of the density n, (A, may be done 
in the LRF and the result is 



Hi (A,£(x)) 



(A) 



(12) 



where rii tCq (A) is the particle density in equilibrium at 
the temperature defined by the parameter A. 
Using Eq. ([12} in ([10} we find 



N; 



and for the ratios 



(A) 



dZ p (x)UK(x) 



(A) 



"3,eq 



(A)' 



(13) 



(14) 
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We thus see that the ratios of hadronic multiplicities 
are exactly the same as the ratios obtained for equilib- 
rium distributions at the temperature A. 

We also note that for the boost-invariant systems we 
have 



Nj IEj-Ay 

dy dy 



(15) 



where Ay is the rapidity range. Hence, in this case, 
the ratios of rapidity densities are the same as the ratios 
obtained for equilibrium at the temperature A. 



IV. TRANSVERSE-MOMENTUM SPECTRA 

In this Section we calculate the transverse-momentum 
spectra in the thermal model with single freeze-out 
and compare our results with the LHC data. We use 
THERMINATOR 2 [H H|[ and choose two options for the 
shape of the freeze-out hypersurface: i) the Cracow 
model, and ii) the blast-wave model. These two mod- 
els have been recently discussed in Ref. [35| . In this 
work we follow closely the method presented in J35| . The 
main difference is the use of the modified distributions at 
freeze-out. 
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A. Cracow model 

In the Cracow model the freeze-out hypersurface 
is defined by the condition that the invariant time 
(t 2 — x 2 — y 2 — z 2 ) 1 / 2 is fixed. The value of this time, 
T3f , and the transverse size of the system, r max , are the 
two geometric parameters of the model. The flow of mat- 
ter at freeze-out has the Hubble form — x^/t^. 

If the hadronic system is in equilibrium, the freeze- 
out conditions arc defined by the two extra parameters: 
the freeze-out temperature, T, and the freeze-out baryon 
chemical potential, fi. In this work we neglect the chemi- 
cal potential and introduce the scale parameter A instead 
of T . In addition, the hadron distributions are taken as 
RSF with the anisotropy parameter £, which we keep 
constant in this Section. Thus, altogether we have four 
independent parameters: T"3f,r max ,A and £. 

In Ref. [35J we analyzed Pb+Pb collisions at the beam 
energy ^/snn = 2.76 TeV for three centrality classes. We 
used three parameters: T3f,r max and T. The freeze-out 
temperature was always set equal to To = 165.6 MeV, 
while the geometric parameters depended on the central- 
ity class, we used: t® { — 9.0 fm and r ^ = 11.4 fm for 
c = 0%-5%, T 3 ° f = 7.4 fm and r° nax = 9.6 fm for c = 10%- 
20%, and finally r 3 ° f = 5.9 fm and r° lax = 7.25 fm for 
c = 30%-40%. 

We have added the superscript to mark that the 
above values have been used in the equilibrium calcu- 
lations. When switching from the equilibrium distribu- 
tions © to the anisotropic distributions ([5]), we use the 



FIG. 1: (Color online) Transverse-momentum spectra of pos- 
itive pions (dots) and kaons (triangles) measured in Pb+Pb 
collisions at the beam energy -^/Inn = 2.76 TeV compared 
to the Cracow model results with £ = 2.5 (a) and £ = —0.5 
(b). The anisotropic hadron distributions are used in the 
model calculations. The equilibrium parameters have been 
rescaled according to the formula (|16[1 . The experimental and 
model results are shown for the centrality class c = 0% — 5%. 

following scheme 

A = Xi 

m = ^(l + O 176 , 

w = ^(1 + o 1/6 . ( 16 ) 

The model transverse-momentum spectra obtained for 
this choice of the parameters are shown in Fig. [T] as solid 
lines and compared to the data. The part (a) describes 
the results for £ = 2.5, while the part (b) shows the 
results for £ = —0.5. The results of the standard ther- 
modynamic fit, obtained in Ref. [351 . are represented by 
the dashed lines. We clearly see that the change intro- 
duced by the momentum anisotropy, induced by the term 
£(p • V) 2 in the exponential function, is very well com- 
pensated by the renormalization of the size of the system 
according to Eq. (TTHl) . 

We note that the power 1/6 appearing in Eq. (fTB")) im- 
plies that the volume of the system in the Cracow model 
scales as (1 + S,) 1 ^ 2 , which compensates the same factor 
in the denominator of Eq. (|12p . Using the same rescaling 
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for T3f and r max we find that the transverse flow profiles 
change very moderately (within 20%). Hence, the invari- 
ance of our results for the spectra with respect to the 
transformation (|16[) is only approximate. This property 
is explained in more detail in the Appendix. 



B. Blast-wave model 

In this Section we present our results obtained with 
the modified version of the blast wave model. This model 
differs from the Cracow model by the form of the freeze- 
out hypersurface and the form of the transverse flow. The 
shape of the freeze-out curve in the Minkowski space is 
controlled by the parameter A, whereas the magnitude 
of the flow is controlled by the parameter Vt- In the 
equilibrium version used in Ref. |35| we have used four 



parameters: A, T , = r" v , and vt- Their values are 
listed in Table 1 of Ref. 

Similarly to the case of the Cracow model, when 
switching from equilibrium to the anisotropic distribu- 
tions, we make the following change of the parameters 



A 

T2i 



To 



a/6 



4(i +0 1 

~° d + n 1/6 

max V 1 S I 



(17) 



Our results are shown in Figs. [5] and [3] Similarly to 
the Cracow model we observe that the spectra are (al- 
most) insensitive to the induced anisotropy provided the 
geometric parameters of the model are appropriately 
rescaled. 



V. HBT RADII 

In Ref. [35[ we have used the Cracow and the blast- 
wave model with different values of the parameter A to 
calculate the pion HBT radii. The main outcome of these 
calculations is that the radii are reproduced best in the 
blast- wave model with A — —0.5. We recall that the 
freeze-out conditions described by a negative value of A 
correspond to the situation where the outer parts of the 
system freeze out earlier than the system's interior. This 
is typical for more advanced hydrodynamic models. 

In this Section we show the results of the calculations 
which are analogous to those presented in Ref. [35j . We 
choose the blast-wave model with A = —0.5 and incor- 
porate the momentum anisotropy at freeze-out. For each 
value of the anisotropy parameter £ we modify the ge- 
ometric parameters according to the formula (JT7J) . As 
we have seen in the previous Section, the modification 
of the geometric parameters compensates the effect of 
the anisotropy in the transverse-momentum spectra. On 
the other hand, the change of the geometric parameters 
should affect the HBT radii, as they are connected with 
the space-time extensions of the system. 
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FIG. 2: (Color online) The transverse-momentum spectra ob- 
tained in the blast-wave model for three different values of the 
parameter A, £ = —0.5 (solid lines) and £ = (dashed lines). 
The equilibrium parameters have been rescaled according to 
the formula (JT7J). The data are taken from |44|. 



Our results for the pion HBT radii obtained in the 
blast- wave model with A = —0.5 for several values of the 
anisotropy parameter £ are shown in Fig. 0] (£ = —0.5 - 
small dots, £ = - solid line, £ = 0.5 - large dots, £ = 1.0 
- long dashes, £ = 2.5 - short dashes). As expected, 
R s ide grows with increasing values of £. This radius has 
a simple geometric interpretation of the transverse size 
of the system, hence, the increase of i? s idc reflects sim- 
ply the growth of r max . We have found that the radius 
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FIG. 3: (Color online) The same as Fig. [2] but for £ = 2.5 
(solid lines) and £ = (dashed lines). The data are taken 
from HI]. 



-Rout behaves in the similar way. On the other hand, the 
radius i?i on g is practically independent of the anisotropy 
parameter £. The latter behavior may be understood as 
a net result of the two effects: a decrease of the longitudi- 
nal homogeneity length due to the decreased longitudinal 
pressure and an increase of this length due to the overall 
increase of the size of the system. Interestingly, the two 
effects almost completely cancel each other. 

The results represented by the solid lines correspond 
to the local equilibrium studied in Ref. [Hj]. They agree 



FIG. 4: (Color online) The pion HBT radii R ou t (a), R sidc (b), 
and i?i ollg obtained in the blast-wave model with the freeze- 
out slope parameter A — —0.5 for five different values of the 
local momentum anisotropy £. The model results are com- 
pared to the ALICE data [451 ] . 

reasonably well with the HBT data. However, Fig. [U 
shows that the agreement with the data may be improved 
if we introduce a non-zero anisotropy. A much better 
agreement with the data is obtained for £ = 1.0 than 
with £ = 0. 

The general conclusion from the HBT calculations pre- 
sented in this Section is that the anisotropy of the mo- 
mentum may be introduced as an additional characteris- 
tics of the system at freeze-out and its non-trivial value 
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may be used to improve the agreement with the data. 



VI. CONCLUSIONS 

In this paper we have considered locally anisotropic 
momentum distributions of hadrons at freeze-out in rel- 
ativistic heavy-ion collisions. We have taken into ac- 
count the local anisotropy between the longitudinal and 
transverse momenta. Our results show that physical ob- 
servables, such as the ratios of hadron abundances or 
the hadronic transverse-momentum spectra, are in prac- 
tice indistinguishable from those obtained in an analo- 
gous equilibrium calculations — the effect of the momen- 
tum anisotropy may be compensated by an appropriate 
change of the geometric models of the system. We have 
also demonstrated that this freedom of the parameters 
may be used to improve the agreement of the model cal- 
culations with the measured HBT radii. 

Our results indicate insensitivity of the thermal ap- 
proach against specific variations of the model assump- 
tions and show that it may be quite successful even in 
the situations where matter is out of equilibrium. 
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VII. APPENDIX 

In this Section, we explain the approximate scaling of 
the transverse-momentum spectra. As an example, we 
consider the Cracow model. In the case of the blast- 



wave model, the arguments leading to the approximate 
scaling are similar. 

In the Cracow model, the spectra of primordial par- 
ticles at zero rapidity are given by the integrals of the 
form 



dN 



dycPpA 



'3f 



(18) 



2 71 



(2tt) s 



dr) 



dd i cosh ■& i sinh ■& i 



x [m^cosh^j^coshry — pj_sinh$j_ cos <fr] 

x exp [— /3(m^cosh?9j_ costly — p±smh , d± cos</>)] . 

Here /3 = l/T, <fi is the azimuthal angle, 77 is the space- 
time rapidity, and 9± is the transverse rapidity connected 
with the transverse distance, r = sinh 6^. 

The momentum anisotropy is implemented by replac- 
ing T by A and by the change of the argument of the 
exponential function, 

m_i_cosh$_i_cosh77 — p_i_sinh!9_L cos <j) — > 



(rn^cosh^coshry — p^sinhz?^ cos < 



£m^_ sinh 77. 



The expression under the square root may be rewritten 
in the equivalent form as 



cos <fi cosh 9± — m±_ cosh 77 sinh 1 



+m 2 + sin 2 4> + Tn]_(l + £) sinh 2 77. 



(19) 



where cosh 77 is relatively suppressed by sinh#^. Since 
the main contribution to the integral (p~8j) comes from 
the region where 77 « 0, the hyperbolic functions may 
be replaced by their approximations, cosh 77 w 1 and 
sinh 77 sa 77. Then, one can easily see that the factor 
1 + £ may be eliminated by the appropriate change of the 
integration variable, which leads to the overall change of 
the normalization of the spectra by the factor (l + £) -1 / 2 . 
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